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Single-spin detection is one of the important challenges facing the develop- 

Q^l ment of several new technologies, e.g. single-spin transistors and solid-state 

Ch ' quantum computation. Magnetic resonance force microscopy with a cyclic adi- 

qh| abatic inversion, which utilizes a cantilever oscillations driven by a single spin, 



is a promising technique to solve this problem. We have studied the quantum 



H ' dynamics of a single spin interacting with a quasiclassical cantilever. It was 
found that in a similar fashion to the Stern-Gerlach interferometer the quantum 
dynamics generates a quantum superposition of two quasiclassical trajectories 
of the cantilever which are related to the two spin projections on the direction 
of the effective magnetic field in the rotating reference frame. Our results show 
that quantum jumps will not prevent a single-spin measurement if the coupling 
between the cantilever vibrations and the spin is small in comparison with the 
amplitude of the radio-frequency external field. 
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Modern solid-state technologies are approaching the level at which manipulating a single 
electron, atom, electron or nuclear spin becomes extremely important. The future success- 
ful development of these technologies depends significantly on development of single-particle 
measurement methods. While a single electron charge can be detected using a single-electron 
transistor, methods for detection of a single electron (or nuclear) spin in solids are still not 
available. However, many proposals for solid-state nano-devices require a single-electron (or 
nucleus) spin measurement. There are a few proposals for a solid-state single-spin measure- 
ment based on the swap operation of a spin state to a charge state [jll, scanning tunneling 
microscopy [Q, or magnetic resonant force microscopy (MRFM) @,|]. 

In this paper, we consider a MRFM single-spin measurement. One of the most promising 
MRFM techniques is based on the cyclic adiabatic inversion (CAI) of electron or nuclear 
spins 0. In this technique, the frequency of the spin inversion is in the resonance with 
the frequency of the mechanical vibrations of the ultrathin cantilever, which allows one to 
amplify the extremely weak force of a spin on the cantilever. The CIA MRFM method was 
successfully implemented as an alternative to the electron and nuclear magnetic resonance 
for macroscopic ensembles of spins It has already achieved a sensitivity, which is 

equivalent to detection of approximately 200 polarized electron spins [0. 

It is clear that the resonant amplification of the cantilever vibrations by a single spin 
cannot be considered as a classical process. Indeed, the driving force acting on the cantilever 
is quantized, since it is proportional to the spin projection. As a result, quantum jumps 
could appear in the cantilever motion, which might prevent the resonant amplification of 
the cantilever oscillations. In fact, the problem of quantum jumps is a very general one. It 
always arises in a continuous observation of a single quantum particle. Despite extensive 
study [^, the quantum jumps in a mechanical motion of classical detectors and their effect 
on a measurement of a quantum system have not been investigated. The analysis of these 
quantum effects and their influence on a single-spin detection in CAI MRFM are the main 
subjects of this paper. 

We consider the cantilever-spin system shown in Fig. 1. A single spin {S = 1/2) is 
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placed on the cantilever tip. The tip can oscillate only in the z-direction. The ferromagnetic 
particle, whose magnetic moment points in the positive z-direction, produces a non-uniform 
magnetic field which acts on the spin. The uniform magnetic field, Bq, oriented in the 
positive z-direction, determines the ground state of the spin. The rotating radio frequency 
(rf) magnetic field, Bi, induces transitions between the ground state and the excited state 
of the spin. The origin is chosen to be the equilibrium position of the cantilever tip with no 
ferromagnetic particle. The r/ magnetic field can be written as Bx = Bi cos[ujt + (f{t)], By = 
—Bism[ujt + f{t)], where (f{t) is a periodic change in phase with the frequency of the 
cantilever, required for a CAl of the spin. A non-uniform magnetic field produces a force 
on the spin which depends on the spin direction. If the spin direction is changed with a 
frequency which equals to the cantilever resonant frequency, the amplitude of the cantilever 
vibrations increases so that it can be detected by optical methods. 

In the reference frame rotating with the frequency of the transversal magnetic field, 
u + dif/dt, the Hamiltonian of the "cantilever-spin" system is, 

p2 ^*yz' / dB 

In Eq. (|l]), Z is the coordinate of the oscillator which describes the dynamics of the quasi- 
classical cantilever tip; Pz is its momentum, m* and Uc are the effective mass and the 
frequency of the cantilever (the mass of the cantilever is: rric = 4m*); Sz and Sx are the z- 
and the x-components of the spin; ul = 'jBz (at z=0) is its Larmor frequency; Ui = 'jBi 
is the Rabi (or nutation) frequency; 7 = gfi/h is the gyromagnetic ratio of the spin; g and 
/i are the g-factor and the nuclear magneton. (We consider for definiteness a nuclear spin, 
but the results can be applied also to an electron spin.) 

We assume that u = ul, which means that the average frequency of the rf field, u, is equal 
to the Larmor frequency of the spin in the permanent magnetic field. Using the dimensionless 
variables: r = uj. and z = Z ^Jm*iJ^Jh the spin-cantilever dynamics is described in the 
rotating frame, by the following Heisenberg operator equation, 

d^z/dr^ + z = 27] Sz , dS/dr = [S x B^jf], (2) 
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where -Be// = (e, 0, —dip/dr + 2rjz) is the dimensionless effective magnetic field, e = uji/ojc 
and 7] = gfi{dBz/dZ)/ {2^Jm^ujpi). Thus, our model includes two dimensionless parameters, 
e and rj. The first one is the dimensionless amplitude of the rf field, and the second one 
describes the interaction of a single spin with a cantilever mechanical vibrations due to a 
non-uniform magnetic field produced by the ferromagnetic particle. Note that the term 
dif/dr in -Be// is caused by the phase modulation of the transversal magnetic field. The 
other part of the z-component of the effective magnetic field, 2riz, describes a nonlinear 
effect - a back reaction of the cantilever vibrations on the spin. 

If the adiabatic conditions (Id'^ip/dr'^l <^ e^) are satisfied and the nonlinear effects are 
small in comparison to the effects of the rf field, the average spin is "captured" by the 
effective magnetic field. More precisely, it precesses around the effective magnetic field 
in such a way that the angle between the directions of the average spin and the effective 
magnetic field approximately remains constant. In CAI MRFM, the value of dip/dr changes 
periodically with the period of the cantilever (27r in our dimensionless variables). As a 
results the effective magnetic field and the average spin change their directions with the 
same period. This leads to a resonant excitation of the cantilever vibrations. 

To test our model, we considered the classical limit of the macroscopic number of spins 
and the classical cantilever. Replacing the operators Sx and Sz in Eq. (||) by the sums of 
operators over all spins in the sample and neglecting the quantum correlation effects, we 
obtain the classical equations of motion for the total average spin and for the cantilever. 
We solved these classical equations numerically for the parameters corresponding to the 
experiment with protons in ammonium nitrate . To estimate the amplitude of stationary 
vibrations of the cantilever within the Hamiltonian approach, we consider the time, r = Qc, 
where Qc is the quality factor of the cantilever. We obtained for the amplitude of stationary 
vibrations of the cantilever Zmax ~ 15 nm, which is close to the experimental value in P|, 
Zmax ~ 16 nm. 

Since a cantilever can be considered as a quasi-classical measuring device, one might 
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expect smooth increase of the cantilever amphtude also in a measurement of a single spin. 
However, a single spin z-component can accept only two values, Sz = ±1/2. Therefore, 
smooth resonant vibrations of the quasi-classical cantilever driven by the continuous oscilla- 
tions of Sz seam to violate the principles of quantum mechanics. From this point of view, one 
should face quantum jumps rather than smooth oscillations of the cantilever. Such jumps 
could prevent the increase of the amplitude of the cantilever vibrations and a single-spin 
detection. To resolve this problem we consider the quantum dynamics of the single spin- 
quasiclassical cantilever system. The dimensionless Schrodinger equation can be written in 
the form. 



where \l/(z,S2,r) is a dimensionless spinor, Ti' = Ti/hujc is the dimensionless Hamiltonian. 
The functions '^i^2{z-,t) are the complex amplitudes to find a spin in the state \sz = ±1/2) 
and a cantilever at the point z at time r. 

To describe the cantilever as a sub-system close to the classical limit, we choose the 
initial wave function of the cantilever in the coherent state Namely, it was taken in the 



form, ^1(2,0) = En=o^n{0)\n), ^2(^,0) = 0, and A„(0) = («"/«!) exp(-|a|V2), where 



|n) is an eigenstate of the oscillator (cantilever) without spin. The initial average values of 
the coordinate and momentum are related to a as (z) = -^(a* ± «), {Pz) = "75 — a). In 
order to correspond the classical limit, we took |ap ^ 1. 

For numerical simulation of the single-spin-cantilever dynamics we used the value of the 
interaction parameter 77 = 0.3. Currently this value is not feasible in experiments with a 
nuclear spin, but can be achieved in experiments with a single electron spin. For instance, for 
the cantilever parameters from and the magnetic field gradient reported in |]^, we obtain 
for a single electron spin rj = 0.8. The phase modulation of the r/ field was taken in the form 
dip/dr = -6000 ± 300r, (if r < 20), and d^p/dr = 1000sin(r - 20), (if r > 20), and the rf 
field amplitude, e = 400. For these parameters the effective magnetic field, 2riz, produced 
by the cantilever vibrations on the spin remains small with respect to the amplitude of the 




(3) 
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rf field. 

The numerical simulations of the quantum dynamics reveal the formation of the asym- 
metric quasi-periodic Schrodinger cat (SC) state for the cantilever. Fig. 2 shows the typical 
probability distributions, P{z,t) — 1^^1(2;, t)|^ -|- \^2{z,t)\^, for nine instants of time. Near 
T = 40, the probability distribution, P{z,t), splits in two asymmetric peaks. After this the 
separation between these two peaks varies periodically in time. The ratio of the peak am- 
plitudes is about 1000 for chosen parameters (the probabilities are shown in the logarithmic 
scale) . 

It is clear that the small peak does not significantly infiuence the average coordinate 
of the cantilever. Fig. 3 shows the average coordinate of the cantilever, {z{t)), and the 
corresponding standard deviation, A(t) = [{z"^) — {^)^]^/^. One can see fast increase of the 
average amplitude of the cantilever vibrations, while the standard deviation still remains 
small. This, in fact, is related to the initial conditions of the spin, which was taken in the 
direction of the z-axis. For instance, if the spin initially points in the x-axis (^i(z,0) = 
^2(2^,0)), our calculations show two large peaks with similar amplitudes. 

The two peaks in the cantilever probability distribution, shown in Fig. 2, indicate two 
possible trajectories of the cantilever. As a result of the consequent measurement of the 
cantilever position the system selects one of the two trajectories. The crucial problem is 
the following: Do the two peaks of the SC state of the cantilever correspond to the definite 
spin states? To answer this question we studied the structure of the wave function of the 
cantilever-spin system. First we have found that both functions, ^1(2;, r) and ^2(-2)T), 
contribute to each peak. Fig. 4 illustrates the probability distributions, 1^1(2;, r)p (red) 
and \^2{z,t)\'^ (blue), for nine instants of time. One can see that these functions have 
maxima at the same values of z. Next, we analyzed the structures of the functions, '^i{z, r) 
and ^2 (-2, t). When two peaks are clearly separated we can represent each of these functions 
as a sum of two terms, corresponding to the "big" and "small" peaks, 

*l,2(^,T)=*?,2(^,T)+*f,2(^,T). (4) 
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We have found that with accuracy up to 1% the ratio, r)/\l/2(-z, t) = 

— \E'2(-2) t) = where the k{t) is a real function independent of z. As a re- 

suh, the total wave function can be represented in the form, 

^{z, s,, t) = ^\z, t)x\s,, t) + ^^{z, r)x^(s„ r), (5) 

where x'^i^z,'^) and x''^{sz,t) are spin wave functions, which are orthogonal to each other. 
Eq. (H) shows that each peak in the probability distribution of the cantilever coordinate 
corresponds to a definite spin wave function. We found that the average spin correspond- 
ing to the big peak (x''|'S'|x'') points in the direction of the vector {e,0, —dLp/dr), whereas 
{x'^\S\x^) points in the opposite direction. Note that up to a small term, 2r]z, the vector 
(e, 0, —dif/dr) is the effective magnetic field acting on the spin (see Eq. (0)). The ratio of 
the integrated probabilities (/ P{z, T)dz) for the small and big peaks (~ 10^^ in Fig. 2) can 
be easily estimated as tan^(B/2), where G is the initial angle between the effective field, 
{e,0, —dip/dr), and the spin direction. Therefore by measuring the cantilever vibrations, 
one finds the spin in a definite state along or opposite to the effective magnetic field. Our 
numerical simulations show that starting with such a new initial condition, i.e. when the av- 
erage spin points along or opposite to the effective field, the probability distribution P{z, r) 
shows again two peaks, as in Fig. 2, but the ratio of the integrated probabilities of these 
peaks is much less (~ 10~^). Thus for chosen parameters, the quantum jumps generated 
by a single spin measurement cannot prevent the amplification of the cantilever vibration 
amplitude, and thus the detection of a single spin. 

So far, the described picture reminds the classical Stern-Gerlach effect in which the 
cantilever measures the spin component along the effective magnetic field. An appearance 
of the second peak, even if the average spin points initially in the direction of the effective 
magnetic field, provides a difference with the Stern-Gerlach effect. The origin of this peak 
is a small deviations from the adiabatic motion of the spin even at large amplitude of the 
effective field, and the back reaction of the cantilever vibrations on the spin. The next 
important question is the following: Is it possible to use CAI MRFM to measure the state of 
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a single spin? We studied numerically the phase of the cantilever vibrations when the initial 
spin points along or opposite the direction of the effective magnetic field. Our computer 
simulations show that the phases of the cantilever vibrations for these two initial conditions 
are significantly different. When the amplitude of the cantilever vibrations increases, the 
phase difference for two initial conditions approaches vr. Thus, the classical phase of the 
cantilever vibrations indicates the state of the spin relatively to the effective magnetic field. 
If the spin is initially in the superposition of these two states, it will acquire one of these 
states in the process of measurement. 

In practical applications it would be very desirable to use CAI MRFM for measurement 
of the initial 2;-component of the spin. For this purpose one should provide the initial 
direction of the effective magnetic field to be a 2;-direction. Then, the initial 2;-component 
of the spin will coincide with its component relatively to the effective magnetic field. In our 
computer simulations presented in Figs. 2-4 we have assumed instantaneous increase of the 
amplitude of the r/ field, at r = 0. It causes the initial angle between the directions of the 
spin and the effective magnetic field, G ~ e/\d(p/dT\ ^ 0.07. To eliminate this initial angle 
we simulated the quantum spin-cantilever dynamics for adiabatic increase of the rf field 
amplitude: e = 20r for r < 20, and e = 400 for r > 20. Dependence for dip/dr was taken 
the same as in Figs. 2-4. The results of these simulations are qualitatively similar to those 
presented in Figs. 2-4, but the integrated probability of the small peak was reduced to its 
residual value ~ 10^^. Neglecting this small probability one can provide the measurement 
of the initial 2;-component of the spin. 

We should also mention that the detection of a single electron spin in an atom can be 



used to determine the state of its nuclear spin [jTO|,|rT|. Such a measurement is possible for 
an atom with a large hyperfine interaction in a high external magnetic field, because the 
electron spin frequency of the atom depends on the state of its nuclear spin. 

In conclusion, we have analyzed the quantum effects in the single-spin measurement 
using cyclic adiabatic inversion (CAI) to drive cantilever vibrations in magnetic resonance 
force microscopy (MRFM). We investigated the quasi-classical cantilever interacting with a 
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single spin using Hamiltonian approach. We have shown that the spin- cantilever dynamics 
generates a Schrodinger-cat (SC) state for the cantilever. The two peaks of the probability 
distribution of the cantilever coordinate correspond approximately to the directions of the 
spin along or opposite to the direction of the effective magnetic field, in the rotating frame. 

In this paper, we did not discuss the intriguing possibility of observing the SC state. 
This requires an estimate of the SC life time, which cannot be derived using our Hamiltonian 
which neglects the environment of the cantilever. Instead, we concentrated on a possibility 
of observing the resonant excitation of the cantilever vibrations, driving by a single spin. We 
demonstrated by a direct computation of the average cantilever position and its standard 
deviation as a function of time that the resonant amplification of the cantilever oscillations 
in indeed possible (for considered region of the system parameters), despite the quantum 
jumps of the single spin. In fact, the standard deviation of the cantilever coordinate becomes 
large only when the angle between the initial spin direction and the effective magnetic 
field approaching 7r/2. In this case the SC state appears with approximately equals peaks. 
However after an observation of the cantilever position the system appears in one of the 
peaks, and the following evolution of the cantilever coordinate shows again the resonant 
amplifications with a very small standard deviation. 

We expect that taking into consideration the interaction with an environment will not 
change our conclusion. Such an interaction will cause the decoherence , which transforms 
the SC state into the statistical mixture. It is clear that this effect as well as the thermal 
vibrations of the cantilever (see, for example cannot prevent an observation of the 

driven oscillations of the cantilever if the corresponding rms amplitude exceeds the rms 
amplitude of the vibrational noise. Another effect of the interaction with the environment 
is the finite quality factor, Qc, of the cantilever, which puts the limit on the increase of the 
cantilever vibrations. The stationary amplitude of the cantilever vibrations can be estimated 
in our Hamiltonian approach by putting r = Qc- 

Finally, we mention two other possible techniques for the cyclic spin inversion in MRFM. 
One of them is the standard Rabi technique. This assumes that in our notation dip/dr = 0, 
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and e = 1, i.e. the Rabi frequency equals to the cantilever frequency. This technique seems 
to be simpler that the CAI MRFM. But the amplitude of the rf field, e, must be much 
greater than the effective field produced by the cantilever on the spin 21] z <^ e = 1. In this 
case, the force acting on the cantilever is very small, and the amplification of the driven 
cantilever vibrations requires a long time, i.e. a large cantilever quality factor. Another 
technique assumes the application of short vr-pulses which periodically change the direction 
of the spin in the time-interval, which is very short in comparison to the cantilever period 
13| . If the time interval between successive pulses equals half of the cantilever period. 



this technique provides a resonant amplification of the cantilever vibrations. Testing this 
technique in MRFM experiments is a challenging problem. 
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FIGURES 




FIG. 1. The cantilever-spin system. Bq is the uniform permanent magnetic field; Bi is the 
rotating magnetic field; M is the magnetic moment of the ferromagnetic particle; S is a single spin. 
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FIG. 2. The probability distribution, P(z), for the cantilever position. The values of parameters 
are: e = 400 and r] = 0.3. The initial conditions are: (^;(0)) = —20, {pz{0)) = (which correspond 
to a = -10^2). The same parameters will be used in all figures below. 
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FIG. 3. Cantilever dynamics, (a) Average coordinate of the cantilever as a function of r and 
(b) its standard deviation A(r) = [(z^(r)) - (z(r))2]V2_ 
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FIG. 4. The probability distributions Pi{z,T) = |^j(z,r)p, i = 1 (red) and i = 2 (blue) 
nine instants in time: Tfc = 92.08 + 0.8A;, A; = 0, 1, 8 
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